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A  coupling-analysis-based  model  to  predict  the  capacity  degradation  of  LiFeP04  batteries  under  multi¬ 
stress  accelerated  conditions  has  been  developed.  In  this  model,  the  joint  effect  on  the  battery  capacity 
degradation  of  any  2  out  of  5  stress  factors,  which  include  ambient  temperature,  end  of  discharge  and 
charge  voltage  (EODV  and  EOCV),  and  discharge  and  charge  rate,  is  studied  through  coupling  validation 
tests.  Coupling  generally  exists  among  these  5  stress  factors,  and  the  coupling  intensity  has  a  certain 
relationship  with  the  stress  levels.  There  is  a  critical  stress  level  at  which  the  coupling  can  be  considered 
negligible,  and  when  the  stress  level  goes  higher,  coupling  aggravates  battery  degradation  exponentially. 
Additionally,  the  study  also  indicates  that  battery  life  shows  stronger  sensitivity  to  discharge  rate  and 
EOCV  than  to  charge  rate  and  EODV.  The  developed  capacity  degradation  model  based  on  the  input  of 
real  operating  conditions  and  coupling  intensity  calibration  achieves  error  less  than  1 5%  when  the  cycling 
goes  into  the  stable  decay  period,  and  the  error  converges  gradually  as  the  cycling  continues. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Battery  degradation  is  generally  predicted  using  one  of  two 
approaches  [1,2].  The  first  one  is  based  on  electrochemical  models 
[3,4]  and  focuses  on  the  micro-mechanism  of  battery  degradation 
on  both  positive  [5,6]  and  negative  electrodes  [7-9],  such  as  active 
particle  loss,  metal  sediment  or  SEI  film  accumulation  [10-13]. 

The  other  approach  was  developed  to  describe  experimental 
phenomena,  like  capacity  fade  and  resistance  increase,  by  fitting 
the  cycling  data  [14-16].  Instead  of  investigating  the  sophisticated 
internal  electrochemical  principles  of  battery  cycling,  the  second 
approach  avoids  the  recognition  of  battery  design  or  manufactur¬ 
ing  parameters,  and  only  external  characteristics,  which  can  be 
measured  easily,  are  required. 

For  the  second  prediction  method,  a  significant  amount  of 
research  has  focused  on  the  effect  of  DOD  (depth  of  discharge), 
discharge  rate,  EODV  or  EOCV  on  the  battery  lifespan,  but  most 
results  have  been  described  roughly  by  experimental  curves  of  the 
capacity  fade  or  increase  in  resistance.  Only  a  few  mathematical 
models  have  been  extracted  and  established.  The  most  widely  used 
external  characteristics  model  represents  the  law  between  ambient 
temperature  and  cycling  capacity  loss  [17,18].  Although  an  Ah- 
throughput  model  [19]  takes  the  coupling  effect  of  temperature 
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and  discharge  rate  into  account,  the  other  models  are  mostly  sin¬ 
gle  stress  accelerated;  moreover,  the  models  do  not  account  for  the 
existence  of  coupling  between  different  stresses. 

Therefore,  these  models  are  not  necessarily  accurate  when 
applied  to  cycle  degradation  prediction  under  complex  real  run¬ 
ning  conditions;  for  instance,  under  real  road  conditions  for  electric 
vehicles. 

In  this  paper,  our  proposed  model  has  several  new  features: 
first,  the  most  important  factors  that  affect  battery  degradation, 

1. e.,  temperature,  EODV,  EOCV,  discharge  rate  and  charge  rate, 
are  all  adopted  as  the  stresses  and  are  used  as  the  input  to  the 
model.  Second,  the  coupling  law  between  these  stresses  is  stud¬ 
ied  and  quantified  by  a  long-term  coupling  validation  test.  Finally, 
a  capacity  degradation  model  based  on  stress  coupling  analysis  is 
established  with  error  evaluation  and  factor  sensitivity  analysis. 

2.  Experiment  design 

The  coupling  validation  test  is  carried  out  as  follows.  In  the  first 
part,  a  single  stress  is  accelerated,  whereas  others  are  set  to  the 
standard  stress  levels  defined  in  Table  1.  The  impact  of  that  single 
factor  on  battery  degradation  is  evaluated  based  on  life  characteris¬ 
tics,  such  as  the  capacity  fade  and  resistance  increase.  In  the  second 
part,  two  factors  are  adjusted  simultaneously  to  accelerated  levels, 
and  the  life  characteristics  of  the  experimental  results  are  com¬ 
pared  to  the  results  in  the  first  part.  If  the  capacity  fade  rate  in  the 
two-factor  accelerated  test  was  similar  to  the  sum  of  the  rates  in 
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Table  1 

Standard  levels  of  5  factors. 


Temperature 

Discharge  rate 

EODV 

Charge  rate 

EOCV 

Standard  levels 

30°C 

1/3  C 

2  V 

1/3  C 

3.65  V 

the  two  corresponding  single-factor  accelerated  tests,  the  coupling 
between  these  two  factors  can  be  ignored,  i.e.,  these  two  factors  are 
almost  decoupled.  If  they  are  quite  different,  it  can  be  determined 
that  there  is  strong,  non-negligible  coupling  between  the  two  fac¬ 
tors,  and  a  further  probe  into  the  relationship  between  the  coupling 
factors  is  required  to  reveal  the  strength  of  the  coupling  and  to 
represent  the  coupling  law  in  the  capacity  degradation  model. 

The  ambient  temperature,  discharge  rate,  EODV,  charge  rate 
and  EOCV  are  selected  as  the  input  stresses.  The  overall  dura¬ 
bility  test  procedures  are  as  follows.  First,  the  single-factor  and 
two-factor  tests  are  performed  to  acquire  the  fade  rate  of  capacity 
and  the  increase  in  resistance,  to  validate  the  existence  of  coupling 
between  the  two  factors  and  to  explore  the  relationship  between 
the  coupling  intensity  and  stress  levels.  Second,  several  groups 
of  single-factor  acceleration  tests  with  different  stress  levels  are 
performed  to  identify  the  parameters  in  single-stress  accelerated 
models.  Third,  using  single-stress  accelerated  models  and  the  cou¬ 
pling  intensity  laws,  the  capacity  degradation  model  is  derived 
taking  into  account  the  coupling  between  factors.  Finally,  the  model 
accuracy  is  verified  experimentally,  and  the  factor  sensitivity  of  this 
model  is  analyzed. 

The  following  criteria  are  determined  for  all  the  cycles  in  the 
test: 

(1)  The  battery  capacity  fade:  A  standard  battery  capacity  test 
method  is  adopted  to  ensure  the  comparability  between  the 
capacity  fade  rates  under  different  cycle  conditions.  The  stan¬ 
dard  battery  capacity  test  has  the  following  parameters:  a 
fixed  room  temperature,  the  standard  charging  method  (see  (3) 
below)  charged  to  full,  rest  for  1  h,  then  1  /3  C  constant  discharge 
to  a  cut-off  voltage  of  2.0  V. 

(2)  The  frequency  of  the  standard  capacity  test:  After  every  20 
cycles  (i.e.,  after  cycle  20, 40, 60, 80, . . .),  the  battery  is  switched 
to  room  temperature  and  gets  a  rest  of  12  h,  and  the  standard 
capacity  and  the  DC  resistance  of  the  battery  are  measured 
under  1/3  C  constant  discharge. 

(3)  Charging  method:  In  each  cycle,  the  charge  rate  and  EOCV  are 
chosen  in  accordance  with  the  experimental  conditions  shown 
in  Tables  3  and  4  using  the  CC-CV  (constant-current,  constant 
voltage)  method.  The  constant  voltage  charging  stage  stops 
when  the  current  drops  to  1  /1 0  of  the  constant  charging  current. 
The  constant  current  and  EOCV  of  the  charging  in  the  standard 
capacity  test  are  1  /3  C  and  3.65  V. 

(4)  Rest  time  between  adjacent  cycles:  5  min  in  general.  However, 
in  the  experimental  results  shown  in  Fig.  4(a),  because  the  cycle 
capacity  declines  rapidly  due  to  harsh  cycling  conditions,  the 


Table  2 

The  parameters  of  the  subject  battery. 


Items 

Parameters  and  units 

Nominal  capacity 

Nominal  voltage 

Weight 

Size 

11  Ah 

3.2  V 
<400  g 

27(+1.5,-0)mm  x  70(±0.2)mm  x  104(±0.4)mm 

Table  3 

Single  factor  accelerated  experiments. 

Test  No. 

Temperature 

Discharge  rate  EODV 

Charge  rate 

EOCV 

1-1 

1-2 

1-3 

1-4 

1-5 

45  °C 

a 

4C 

1.25  V 

1.5  C 

3.95  V 

a  The  blanks  refer  to  ‘standard  stress  levels’  from  Table  1 ;  the  same  system  is  used 
below. 


rest  time  is  changed  from  5  min  to  2  h,  and  the  adjustment  point 
is  shown  in  the  figure. 

(5)  When  to  stop  cycling:  Cycling  continues  until  one  of  the  fol¬ 
lowing  occurs:  the  standard  battery  capacity  drops  to  60%  of 
the  nominal  capacity,  the  cycling  capacity  drops  to  1 5%  of  the 
nominal  capacity,  a  total  of  200  cycles  is  completed,  or  unsafe 
circumstances  occur,  such  as  battery  bulging  or  a  sharp  rise  in 
temperature. 

(6)  Temperature  adjustment:  After  the  adjustment  of  the  ambient 
temperature,  the  battery  should  first  rest  for  more  than  12  h 
until  the  internal  and  external  temperatures  coincide. 

(7)  Resistance  test  method:  HPPC  test  [20]  with  a  pulse  current 
of  1  C. 

2. 1 .  Single  factor  accelerated  experiments 

The  experiments  were  performed  on  some  new 
LiFeP04/graphite  cells  and  the  parameters  of  the  battery  are 
shown  in  Table  2.  For  all  the  experiments  in  Table  3,  only  one  stress 
is  accelerated,  and  the  others  remain  at  the  standard  levels  shown 
in  Table  1. 

2.2.  Double-factor  accelerated  experiments 

The  experiments  were  carried  out  on  1 0  cells  of  the  same  kind  as 
those  in  Section  2.1 .  Each  of  the  experiments  has  double  accelerated 


Table  4 

Double-factor  accelerated  experiments. 


Test  No. 

Temperature 

Discharge  rate 

EODV 

Charge  rate 

EOCV 

2-1 

45  °C 

4C 

2-2 

45  °C 

1.25V 

2-3 

45  °C 

1.5  C 

2-4 

45  °C 

3.95  V 

2-5 

4C 

1.25V 

2-6 

4C 

1.5  C 

2-7 

4C 

3.95  V 

2-8 

1.25V 

1.5  C 

2-9 

1.25V 

3.95  V 

2-10 

1.5  C 

3.95  V 
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Table  5 

The  calculation  of  coupling  error. 


Single  stress 

Corresponding  double  stress 

Coupling  error  A  (%) 

Test  No. 

SCFR 

Test  No. 

SCFR 

Test  No. 

SCFRd 

1-1 

0.00576 

1-2 

0.01425 

2-1 

0.02025 

-1.19 

1-1 

0.00576 

1-3 

0.00543 

2-2 

0.00527 

112.33 

1-1 

0.00576 

1-4 

0.00474 

2-3 

0.00798 

31.58 

1-1 

0.00576 

1-5 

0.00481 

2-4 

0.00535 

97.57 

1-2 

0.01425 

1-3 

0.00543 

2-5 

0.02221 

-11.39 

1-2 

0.01425 

1-4 

0.00474 

2-6 

0.02065 

-8.04 

1-2 

0.01425 

1-5 

0.00481 

2-7 

0.01783 

6.90 

1-3 

0.00543 

1-4 

0.00474 

2-8 

0.00824 

23.42 

1-3 

0.00543 

1-5 

0.00481 

2-9 

0.00826 

23.97 

1-4 

0.00474 

1-5 

0.00481 

2-10 

0.00705 

35.46 

stresses,  and  these  10  sets  in  Table  4  cover  every  combination  of  2 
stresses. 

The  stress  levels  in  test  2-1  to  2-10  correspond  to  those  in  test 
1-1  to  1-5.  For  example,  the  two  accelerated  stress  levels  in  2-1 
are  designed  to  be  the  combination  of  accelerated  stress  levels  in 
1-1  and  1-2  to  explore  the  relationship  between  temperature  and 
discharge  rate. 


2.3.  Parameter  recognition  experiments  of  single -factor  models 

Further  experiments  were  carried  out  on  6  new  cells  of  the 
same  kind  as  those  in  Section  2.1  to  identify  the  parameters  in  Eqs. 
(3)-(7). 


3.  Experiment  results 

3.1.  Results  of  the  single  factor  accelerated  experiments 


x  =  0.  Thus,  the  standard  capacity  fade  rate  (SCFR)  is  defined  as 
follows:  /=  (a  -y)/x,  with  a  unit  of  Ah( cycle)-1 . 

The  SCFR  of  all  single-factor  experiments  are  shown  in  Table  5. 

3.2.  Results  of  the  double-factor  accelerated  experiments 

The  cycling  curves  of  test  2-3  in  Table  4  are  shown  in  Fig.  3 
to  represent  all  the  double-factor  accelerated  experiments.  Only 
partial  curves  are  plotted  due  to  the  large  number  of  cycles. 

The  cycling  capacity,  the  standard  capacity  and  the  resistance 
for  each  test  in  Table  4  are  shown  in  Fig.  4. 

The  SCFR  of  all  double-factor  experiments  are  shown  in  Table  5. 

Test  2-1  in  Table  4  only  cycled  60  times  due  to  a  sharp  decline 
in  its  cycling  capacity  under  45 °C  with  a  4C  discharge  rate.  After 
extending  its  rest  time  from  5  min  to  2  h,  the  rapid  decay  eased,  as 
shown  in  Fig.  4(a).  However,  after  60  cycles,  the  cycling  capacity 
remained  less  than  2  Ah.  To  ensure  the  safety  of  laboratory  person¬ 
nel  and  equipment,  the  cycling  experiment  was  stopped  after  the 
60th  cycle. 


The  cycling  curves  of  test  1-2  in  Table  3  are  shown  in  Fig.  1  as 
an  example  of  all  the  single-factor  accelerated  experiments.  Only 
parts  of  curves  are  plotted  because  many  cycles  were  completed. 

The  capacity  of  each  cycle  and  the  standard  capacity  after  every 
20  cycles  can  be  calculated  according  to  the  cycling  discharge 
curves  and  the  standard  capacity  tests,  while  the  ohmic  resistance 
and  total  resistance  can  be  measured  by  the  HPPC  method  after 
every  20  cycles.  The  cycling  capacity,  the  standard  capacity  and 
the  resistances  are  shown  in  Fig.  2,  and  the  ohmic  resistance  and 
total  resistance  during  discharge  and  charge  are  noted  as  Rdch-ohm> 
Rdch-tot.  Rcha-ohm  and  Rcha-tot  in  the  figure. 

The  fitted  line  of  each  standard  capacity  curve  in  Fig.  2  is 
assumed  to  have  the  form  y  =  a-fx.  This  equation  describes  the 
relationship  between  the  standard  capacity  y,  the  cycling  times  x 
and  the  capacity  fade  rate  f  where  a  is  the  initial  capacity  when 


4.  Analysis  of  the  experimental  results 


The  slopes  of  linear  fit  lines  to  the  standard  capacity  degradation, 
namely,  the  SCFR,  are  listed  in  Table  5  with  a  unit  of  Ah(cycle)-1 . 
The  coupling  errors  in  Table  5  are  calculated  by  Eq.  (1): 


E l,SCFRi-SCFRd 
SCFRd 


x  100% 


(1) 


where  J^=1SCFRj  t^ie  sum  °f  the  SCFR  values  of  2  single-factor 
tests  and  SCFRd  is  the  SCFR  of  the  corresponding  double-factor  test 
listed  in  the  same  row  in  Table  5. 

The  coupling  errors  in  Table  5  reveal  significant  information 
about  the  joint  effect  of  these  5  influential  factors: 


Time/  s 

(a)  Cycling  discharge  curves 


(b)  Standard  capacity  tests 


Fig.  1.  Test  No.  1-2  cycling  discharge  curves  (partial)  and  standard  capacity  tests. 
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(a)  Test  No.  1  - 1  (b)  Test  No.  1  -2 


Cycle 

(c)  Test  No.  1-3 


Cycle 

(d)  Test  No.  1-4 


3 


Cycle 

(e)  Test  No.  1-5 

Fig.  2.  (a)-(e)  The  cycling  capacity,  the  standard  capacity  and  the  resistances  of  single  factor  accelerated  experiments. 


Time/  s 


(a)  Cycling  discharge  curves 


(b)  Standard  capacity  tests 


Fig.  3.  Test  No.  2-3  cycling  discharge  curves  (partial)  and  standard  capacity  tests. 


Resistance/  mfi  Resistance/  mfl 
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Cycle  Cycle 


(a)  Test  No.  2-1 


(b)  Test  No.  2-2 


(c)  Test  No.  2-3 


Cycle 

(d)  Test  No.  2-4 


(e)  Test  No.  2-5  (f)  Test  No.  2-6 


60 

50 

40  m 
w 

30  % 
8 

20| 

10 

0 


Cycle 

(g)  Test  No.  2-7 


Cycle 

(h)  Test  No.  2-8 


©Test No.  2-10 


Fig.  4.  (a)-(j)  The  cycling  capacity,  the  standard  capacity  and  the  resistance  of  double-factor  accelerated  experiments. 


Resistance/  mQ 
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(a)  The  reverse  trend  (b)  The  exponential  fit 

Fig.  5.  The  relationship  between  A  and  SCFRd. 


(1)  The  value  of  the  coupling  error  reflects  the  error  of  the 
entirely-decoupled  assumption  compared  with  the  real  situ¬ 
ation.  Table  5  shows  that  there  is  coupling  between  almost 
every  set  of  two  stresses.  Neglecting  coupling  considerably  will 
reduce  the  model  accuracy  due  to  the  relatively  large  coupling 
errors,  as  shown  in  the  table. 

(2)  The  A  and  SCFRd  in  Table  5  demonstrate  an  obvious  reverse 
trend  and  fit  an  exponential  decay  relationship  with  a  coef¬ 
ficient  of  determination  R2  =  0.94,  as  shown  in  Fig.  5(a)  and 
(b).  Although  the  determination  coefficient  of  the  fitting  is  not 
perfect,  because  it  is  impossible  (and  not  required)  to  forecast 
the  battery  lifespan  as  accurately  as  1  month  or  several  days, 
the  accuracy  of  the  battery  degradation  prediction  model  is  far 
below  the  general  requirements  in  dynamic  controls.  Therefore, 
this  fitting  coefficient  is  acceptable. 

Based  on  the  lack  of  one-time  damages,  as  SCFRd  increases, 
the  coupling  error  approaches  a  negative  value.  The  coupling 
error  equals  zero  at  a  critical  SCFRd ;  thus,  the  two  factors 
are  completely  decoupled  and  affect  battery  degradation  sep¬ 
arately.  When  the  SCFRd  is  less  than  the  critical  value,  the 
coupling  error  is  positive,  and  the  interaction  between  the  two 
factors  reduces  the  degradation  of  the  battery.  When  the  SCFRd 
is  greater  than  the  critical  value,  the  coupling  error  is  nega¬ 
tive,  and  the  degradation  becomes  more  intense  due  to  the 
interaction  between  these  two  factors.  Therefore,  the  impact 
of  coupling  on  battery  capacity  can  be  obtained  from  the  stress 
level  input. 

(3)  When  the  stress  levels  are  beyond  the  normal  scope,  the  capac¬ 
ity  fading  mechanism  might  change;  thus,  Fig.  5(b)  should  be 
applied  in  a  limited  range.  The  applicable  range  is  derived  in 
the  following  modeling  process. 

5.  Model  development 

5.1.  Model  derivation 

For  a  battery  capacity  degradation  model  with  five  independent 
stresses  as  input,  the  SCFR  under  joint  effect  can  be  expressed  as 
follows: 


where  SCFRT,  SCFRdr,  SCFReodv,  SCFRCr  and  SCFRE0Cv  refer  to  the 
SCFR  when  the  single  stress  accelerated  input  is  temperature,  dis¬ 
charge  rate,  EODV,  charge  rate  and  EOCV,  respectively. 

Several  single  stress  durability  models,  such  as  the  Arrhenius 
and  Inverse  Power  Law  models  [21,22],  have  been  developed  and 
can  be  applied  to  describe  the  relationship  between  battery  lifespan 
and  single  factors: 


h 

SCFRt  =  — 

ST 


b 

ATecr/T 


(3) 


SCFRdr  =  A 

SDR 


b 

BdrIdr-01* 


(4) 


SCFREqdv  = 


b 

%EODV 


b 

Beodv  vEODV  ~°eodv 


(5) 


SCFRcr  =  A 

SCR 


b 

Bcr/cr"0cr 


(6) 


E0CV  %EOCV  B EOCV  V EOCV  Deocv  ^  ^ 

where  i;T,$DR,  %eodv>  §cr»  and  %eocv  are  the  cycle  life  when  the  sin¬ 
gle  stress  accelerated  input  is  temperature,  discharge  rate,  EODV, 
charge  rate  and  EOCV,  respectively;  AT,  BDR,  BEOdv ,  Bcr ,  and  BEOcv 
are  the  constant  parameters  before  the  power  function,  and  CT,  Ddr, 
DEodv ,  Dcr>  and  DEOcv  are  the  constant  parameters  in  the  power 
function.  T is  the  temperature  in  K;  IDR  and  ICR  are  the  discharge  and 
charge  rate;  vEOdv  and  vEOcv  are  the  EODV  and  EOCV;  b  is  defined  as 
b  =  a  —  y sec y  in  which  ysec  is  the  standard  capacity  when  the  battery 
life  ends. 

Thus,  Eq.  (2)  can  be  substituted  by: 


b  (  1  1  1 _ 

SCFRq  \AEecT/T  BdrIdr~Udr  BEodvVeodv~Deodv 


1 

+  BCr/Cr-°cr 


+ 


_ i _ y1 

B EOCV  V EOCV  Deocv  ) 


(3) 


where  §0  is  the  cycle  life  under  five  independent  accelerated 
stresses. 

Then,  this  model  is  modified  to  take  into  account  the  coupling 
between  stresses  according  to  the  coupling  error  analysis  per¬ 
formed  in  Section  4: 


A  A 

SCFR0c  =  -f  (SCFRT  +  SCFRdr)  +  -f(SCFRT  +  SCFRE0DV)  +  -f(SCFRT  +  SCFRcr)  +  -f[SCFRT  +  SCFRE0Cv) 

+  d-f  (SCFRdr  +  SCFREOdv)  +  SCFRdr  +  SCFRcr)  +  SCFRDr  +  SCFReocv)  +  SCFRE0DV  +  SCFRcr)  (9) 

+^SCFReodv  +  SCFReocv )  +  A(SCFRCr  +  SCFRE0CV) 


SCFRq  =  SCFRj  +  SCFRdr  +  SCFREqdv  +  SCFRcr  +  SCFREqcv 


where  SCFR0C  is  the  SCFR0  after  considering  stress  coupling.  3i, 
(2)  82  . . .  9io  are  the  coupling  correction  factors  between  every  2 
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Fig.  6.  The  relationship  between  coupling  error  and  the  two  single-factor  SCFR. 


stresses.  Each  of  the  factors  equals  1  and  Eq.  (9)  is  transformed 
to  Eq.  (2)  under  the  assumption  of  complete  decoupling. 

An  example  of  how  to  calculate  3i  is  given  to  illustrate  the  cal¬ 
culation  method  for  all  the  coupling  correction  factors. 


(1 )  is  the  coupling  correction  factor  of  the  stress  of  temperature 
and  discharge  rate,  and  Eq.  (1 )  is  applied  to  obtain  the  coupling 
error  of  this  combination: 


Aai 


SCFRt  +  SCFRdr  -  SCFRt-DR 
SCFRj-dr 


x  100% 


(10) 


where  SCFRT-dr  is  the  SCFRd  with  temperature  and  discharge 
rate  as  double  accelerated  stresses.  SCFRT  and  SCFRdr  in  this 
formula  is  calculated  from  Eqs.  (3)  and  (4). 

(2)  Eq.  (11),  which  is  obtained  from  the  exponential  fitting  of 
Fig.  5(b),  provides  another  formula  for  Aai  and  SCFRT-dr  in 
addition  to  Eq.  (10). 

SCFRT-DR  =  0.01306e-A9i/°'28794  +  0.00  4  22  (11) 


The  relationship  between  Aai  and  {SCFRt  +  SCFRdr)  derived 
from  Eqs.  (10)  and  (11)  is  shown  in  Fig.  6. 

In  Fig.  6,  when  A  is  in  the  range  of  [-0.7,0.8],  A 
and  {SCFRt  +  SCFRdr)  have  one-to-one  correspondence,  and 
{SCFRt  +  SCFRdr)  is  in  the  range  of  [0.0091, 0.0458]. 

When  {SCFRt  +  SCFRdr  )  e  [0.0091,  0.0173) 

=*  Ae(0,  0.8] 

=K SCFRt  +  SCFRdr  )  >  SCFRT-dr ,  the  coupling  between  factors 
lessens  the  battery  capacity  fade. 

When  {SCFRt  +  SCFRdr)  e  (0.0173,  0.0458] 

=>  A  e  [-0.7,0) 

=^{SCFRt  +  SCFRdr)  <  SCFRT-dr ,  the  coupling  between  factors 
speeds  up  the  battery  capacity  fade. 

When  the  stress  is  too  large  or  small,  the  sum  of  two  single¬ 
factor  SCFR  may  fall  outside  the  one-to-one  correspondence 
range.  If  the  stress  is  not  far  off  the  one-to-one  interval,  its 
corresponding  coupling  error  can  be  calculated  through  linear 
extrapolation  based  on  the  start  and  end  points  of  the  interval.  If 
the  deviation  is  higher,  the  battery  fading  mechanism  may  have 
been  changed,  and  the  above  rules  may  be  no  longer  applicable. 

Because  the  range  outside  the  one-to-one  interval  exceeds 
the  experimental  scope  in  Section  3,  it  can  only  be  concluded 
that  if  the  stress  is  in  the  right  lateral  region  of  the  one-to- 
one  interval,  the  stress  level  is  extremely  small,  the  coupling 
between  factors  significantly  reduces  the  battery  fade,  and  the 
SCFR  under  double-factor  conditions  is  much  smaller  than  the 
sum  of  the  fading  rates  of  the  two  single-factor  experiments. 
Therefore,  the  coupling  error  is  so  big  that  the  correction  factor 
in  Eq.  (12)  is  close  to  zero.  Similarly,  if  the  stress  is  in  the  left 


Table  6 

Parameter  recognition  results  for  single-factor  Eqs.  (3)-(7). 


SCFR 

Multiplier  factor  parameters 

Power  parameters 

SCFRt 

Ar  =  0.0635 

Cr  =  2838.0121 

SCFRdr 

Bdr  =  624.8817 

Ddr  =  0.6905 

SCFReodv 

Beodv  =  547.8087 

Deodv  =  0.8514 

SCFRcr 

Bcr  =  641.6489 

Dcr  =  0.4754 

SCFReocv 

#£0cv  =  7.6638  x  105 

Deocv  =  5.1960 

Table  7 

Model  accuracy  test. 

Test  parameter  items 

Parameter  value  and  unit 

Ambient  temperature 

50  °C 

Discharge  rate 

1/3  C 

EODV 

2  V 

Charge  rate 

1/3  C 

EOCV 

4.2  V 

Cycling  times 

110 

SFCR 

0.01337 

lateral  region,  the  coupling  error  tends  to  - 1 ,  and  the  correction 
factor  tends  to  positive  infinity. 

Given  the  stress  levels  of  each  of  two  single  factors,  the  cou¬ 
pling  error  can  be  calculated  by  the  look-up  table  provided  in 
Fig.  6. 

(3)  Since  9i  is  defined  as  the  ratio  of  SCFR  under  double-factor 
conditions  to  its  corresponding  sum  under  two  single-factor 
conditions,  it  can  be  written  as  follows: 

_  SCFRt-dr  _  1  /px 

1  “  SCFRt  +  SCFRdr  ~  1  +  Aai  1  J 

(4)  9i  can  be  solved  using  Eqs.  (10),  (11)  and  (12)  together.  Similar 
methodology  can  be  applied  to  the  derivation  of  other  cou¬ 
pling  correction  factors,  and  SCFR £  in  Eq.  (9)  can  be  calculated  to 
predict  the  cycle  life  under  multi-stress  accelerated  conditions: 


*s  = 


b 

SCFRq 


(13) 


5.2.  Parameter  recognition  and  model  accuracy  validation 


Model  parameters  are  identified  and  the  model  accuracy  is  ver¬ 
ified  through  cycling  tests  in  this  section. 

First,  the  battery  life  is  assumed  to  end  when  the  standard  capac¬ 
ity  declines  to  70%  of  nominal  capacity.  Thus, 

b  =  a-ysec  =  11  x  (1  -  70%)  =  3.3  Ah 

Second,  several  extra  single-factor  experiments  are  conducted 
to  determine  the  SCFR,  which  can  be  used  in  the  model  parame¬ 
ters  identification  of  Eqs.  (3)-(7).  The  parameters  in  Eqs.  (3)-(7) 
determined  from  these  experiments  are  shown  in  Table  6. 

Then,  a  test  is  carried  out  with  the  conditions  shown  in  Table  7 
to  verify  the  model  accuracy. 

If  each  stress  in  Table  7  occurs  while  the  other  stresses  remain 
at  the  standard  levels,  the  SCFR  of  these  five  single-factor  acceler¬ 
ated  experiments  are  counted  in  Table  8  according  to  the  parameter 
recognition  results  in  Table  6  and  the  inputs  in  Table  7. 

The  corresponding  coupling  errors  and  correction  factors 
between  the  factors  in  Table  7  are  calculated  and  shown  in  Table  9 
using  the  look-up  curve  in  Fig.  6. 


Table  8 

SCFR  of  single-factor  accelerated  experiments  under  the  input  of  Table  7. 


SCFRt 

SCFRdr 

SCFRdr 

SCFRcr 

SCFReocv 

0.0079 

0.0025 

0.0033 

0.0030 

0.0075 
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Table  9 

The  coupling  errors  and  correction  factors  using  the  input  values  shown  in  Table  7. 


Correction 
factor  symbol 

Sum  of  two 
single-factor  SCFR 

Coupling 

error 

Correction 

factors 

9i 

0.0104 

0.41 

0.7092 

d2 

0.0112 

0.32 

0.7576 

83 

0.0109 

0.35 

0.7407 

d4 

0.0154 

0.07 

0.9346 

d5 

0.0058 

NAa 

0 

d6 

0.0055 

NA 

0 

d? 

0.01 

0.47 

0.6803 

d8 

0.0063 

NA 

0 

d9 

0.0108 

0.36 

0.7353 

8 10 

0.0105 

0.4 

0.7143 

a  Both  of  the  stress  levels  are  quite  small,  and  the  stress  level  is  to  the  right  of  the 
one-to-one  interval;  thus,  the  correction  factor  can  be  assumed  to  be  approximately 
zero. 


Finally,  the  results  in  Tables  8  and  9  can  be  substituted  into  Eq. 
(9),  and  the  SCFR  under  the  joint  effect  of  5  factors,  which  has  taken 
into  stress  coupling  into  account,  can  be  calculated  as  follows: 

SCFRcou  =  0.0151. 

The  SCFR  without  considering  stress  coupling  is  counted  by  the 
sum  of  all  single-factor  SCFR  in  Table  8: 

SCFRdecou  =  0.0242. 

Compared  with  the  real  experimental  result  of  the  SCFR  in 
Table  7  ( SCFRreai  =  0.01337),  the  SCFR  value  without  accounting  for 
stress  coupling  ( SCFRdecou )  caused  a  forecast  error  of  81%,  whereas 
the  SCFRCou,  derived  from  the  capacity  degradation  model  that 
has  been  modified  by  stress  coupling,  achieved  an  error  less  than 
12.94%.  The  three  SCFR  are  compared  in  Fig.  7. 

With  the  ongoing  cycling  of  the  experiment  in  Table  7,  the 
SCFRreai  obtained  by  the  known  data  is  gradually  changing.  Fig.  8 
shows  that  the  model  error  fluctuates  with  an  increasing  number 
of  cycles.  It  can  be  concluded  that  for  the  first  few  dozen  cycles, 
the  model  error  can  be  up  to  50%  due  to  fluctuation  of  the  capacity 
fading  rate;  when  the  cycle  number  accumulates  to  60  or  more,  the 
absolute  value  of  the  model  error  is  reduced  to  less  than  20%.  When 
the  cycle  number  increases  to  1 1 0,  the  model  error  is  approximately 
13%.  The  error  converges  non-monotonically  as  the  cycling  contin¬ 
ues,  which  is  relevant  to  the  stabilization  of  the  standard  capacity 
fading  rate  as  the  number  of  cycles  increases. 

5.3.  Factor  sensitivity  analysis 

The  capacity  degradation  model  provides  an  analysis  tool  to 
explore  the  sensitivity  of  battery  durability  to  these  five  factors,  and 


Cycle 

Fig.  8.  Error  convergence  of  the  durability  model. 


the  relationship  between  SCFR  and  the  five  factors  can  be  obtained 
from  the  derivative  of  Eq.  (9). 

dSCFRcQ  dSCFRcQ  dSCFRT 
dT  ~  dSCFRT  '  dT 

(di  T  d2  +  d3  +  d4)  b_  ct/t  Ct  r^A\ 

4  'At  T 2  [  J 


dSCFRc0  dSCFRc0  dSCFRDR 
dim  ~  dSCFRDR  '  dIDR 

=  .  Jl.Ddr{Idr)ddr- i 

4  Bdr 


(15) 


3SCFRC0 

dVEODV 


dSCFR q  dSCFREOpy  _  {d2  +  ^5  +  ^8  +  ^9) 
SSCFReodv  dvEODV  4 


b 

Beodv 


(- DeodvXveodv )  Deodv 


(16) 


3SCFRC0  dSCFRc0  dSCFRcR 
dICR  ~  dSCFRCR  '  dICR 

=  (93  +  96+49s  +  9i°)  ■  A  ■  DCR(ICRf™' 


(17) 


dSCFR q  SSCFRq  dSCFRpocv  (94  +  ^7  +  dg  +  9io) 

Bveocv  dSCFRRocv  Bveocv  4 

•  n - •  deocv(veocv)Deocv 

°EOCV 


(18) 


The  following  condition  is  adopted  as  an  example  to  analyze  the 
factor  sensitivity  of  the  durability  model:  ambient  temperature  of 
40  °C,  a  discharge  rate  of  3  C,  EODV  of  1 .5  V,  a  charge  rate  of  1  C,  and 
EOCV  of  4  V.  The  coupling  errors  and  the  correction  factors  under 
this  condition  are  calculated  and  shown  in  Table  10. 

Given  the  parameter  identification  in  Table  6  and  the  correction 
factors  in  Table  1 0,  the  result  of  the  capacity  fading  rate  of  this  sensi¬ 
tivity  analysis  case  is  around  0.0338,  based  on  Eq.  (9).  The  derivative 
results  in  Table  11  are  obtained  by  calculating  Eqs.  (14)-(18). 

Sensitivity  analysis  can  help  determine  the  maximum-weighted 
influence  factors  for  battery  life.  Flere,  the  concept  of  ‘equivalent 
stress  increment’  is  introduced  to  describe  a  group  of  factor  incre¬ 
ments  that  affect  battery  life  equivalently.  The  variables  AT,  A IDR, 
Aveodv,  AIDc  and  AvEqcv  in  the  following  formula  form  a  group  of 
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Table  10 

The  coupling  errors  and  correction  factors  for  the  case  of  factor  sensitivity  analysis. 


Correction 
factor  symbol 

Sum  of  two 
single-factor  SCFR 

Coupling 

error 

Correction 

factors 

0.0173 

0 

1 

d2 

0.0145 

0.12 

0.8929 

0.0111 

0.33 

0.7519 

d4 

0.0118 

0.28 

0.7813 

a5 

0.0198 

-0.08 

1.0870 

de 

0.0164 

0.03 

0.9709 

d7 

0.0171 

0.01 

0.9901 

d8 

0.0137 

0.16 

0.8621 

d9 

0.0143 

0.13 

0.8850 

9io 

0.0109 

0.35 

0.7407 

Table  11 

The  derivative  results  for  factor  sensitivity  analysis. 

dSCFRC 

dSCFR^ 

3 SCFRC 

3 SCFRC 

3  T 

81  DR 

8vEODV 

dIDC 

dvEocv 

1.488  x  10“* 1 2 3 4 

0.0026 

-0.0022 

0.0020 

0.0064 

equivalent  stress  increments  to  battery  life. 

dSCFRcn  dSCFRi  dSCFRcn 

°- AT  =  -^A/D]?  =  3 - °-AvE0Dv 


dT 


dlDR 


dSCFRcn  A 

ICR  = 


dVEODV 

dSCFRcn 


dvEOCV 


-A  vEocv 


(19) 


A  group  of  equivalent  stress  increments,  as  in  Eq.  (20),  are  cal¬ 
culated  from  the  derivatives  in  Table  1 1 . 

r  AT  =  10  °C 

AIe)R  =  0.5722  C 

<  Aveodv  =  -0.6762  V  (20) 

AICr  =  0.7439  C 
<  A Veocv  —  0.2325  V 

The  group  of  equivalent  stress  increments  indicates  that  each 
one  of  the  following  events  will  have  an  equivalent  impact  on  bat¬ 
tery  life:  a  temperature  increase  of  10°C,  a  discharge  rate  increase 
of  0.57  C,  an  EODV  decrease  of  0.68  V,  a  charge  rate  increase  of  0.74  C 
and  an  EOCV  increase  of  0.23  V.  Hence,  the  following  conclusions 
on  the  model  sensitivity  can  be  drawn: 

(1 )  The  ambient  temperature  significantly  affects  battery  life.  Tak¬ 
ing  the  cycling  conditions  of  the  sensitivity  analysis  as  the 
starting  point,  a  10°C  rise  in  temperature  will  increase  the 
capacity  fading  rate  by  0.0015  and  shorten  the  calendar  life  by 
5%.  As  with  cycling  under  high  temperature,  cycling  under  low 
temperatures  also  accelerates  the  degradation  of  batteries.  The 
natural  temperature  range  is  as  wide  as  40  °C  in  Beijing  from 
summer  to  winter,  and  a  proper  heating  and  cooling  system  for 
battery  packs  will  help  achieve  a  longer  battery  life. 

(2)  The  battery  life  will  be  shortened  by  a  too  low  EODV  or  too  high 
EOCV,  and  the  sensitivity  to  EOCV  is  significantly  higher  than 
to  EODV.  In  the  above  example,  the  0.23  V  increase  of  EOCV  has 
an  equivalent  impact  on  battery  life  to  a  0.68  V  drop  of  EODV. 
Compared  to  over-discharge,  over-charge  is  more  likely  to  be  a 
fatal  mistake  for  lithium-ion  batteries,  and  the  charging  process 
should  be  carefully  controlled. 

(3)  Excessive  charge  and  discharge  rates  harm  batteries,  and  bat¬ 

tery  life  is  more  sensitive  to  discharge  rate  than  to  charge  rate. 
In  the  above  example,  increasing  the  discharge  rate  by  0.57  C  is 
equivalent  to  increasing  the  charge  current  by  0.74  C.  Addition¬ 
ally,  the  standard  charge  is  more  favorable  for  batteries  than  the 


fast  charge.  Given  the  same  operating  conditions,  including  the 
same  EOCV,  if  the  standard  charge  rate  (e.g.,  1/3  C)  is  changed 
to  a  fast  charge  rate  (e.g.,  3  C),  the  battery  life  will  be  shortened 
by  more  than  20%. 

It  should  be  noted  that  the  stress  increments  in  formula  (20)  are 
based  on  the  starting  point  of  the  conditions  in  this  sensitivity  anal¬ 
ysis.  The  data  in  Eq.  (20)  and  in  the  above  conclusions  are  subject 
to  change  when  the  starting  point  has  been  altered. 

6.  Conclusions 

A  battery  capacity  degradation  model  with  the  multi-factor 
input  is  studied  in  this  paper.  First,  durability  tests  are  performed  to 
obtain  the  capacity  fade  and  resistance  increase  with  cycling  and  to 
verify  the  general  existence  of  coupling  between  stresses.  Accord¬ 
ing  to  the  coupling  verification  test  results,  the  relevance  between 
coupling  intensity  and  stress  levels  is  further  investigated.  Then, 
the  multi-factor  capacity  degradation  model  is  established  based 
on  the  coupling  intensity  calibration  and  the  input  of  real  operat¬ 
ing  conditions.  This  model  achieves  an  error  less  than  1 5%  when  the 
cycling  goes  into  the  stable  decay  period,  and  the  error  converges 
gradually  as  the  cycling  continues.  Finally,  model  sensitivity  analy¬ 
sis  was  performed  to  determine  the  different  weights  of  these  five 
factors  with  respect  to  battery  life. 

The  main  conclusions  of  this  paper  are  as  follows: 

( 1 )  The  5  influencing  factors  of  battery  capacity  degradation,  which 
are  ambient  temperature,  discharge  rate,  EODV,  charge  rate, 
and  EOCV,  are  coupled  in  general. 

(2)  The  coupling  is  related  to  stress  levels.  There  exists  a  critical 
stress  level  at  which  the  coupling  can  be  relatively  neglected. 
When  the  stress  level  exceeds  the  critical  value,  the  coupling 
aggravates  the  battery  degradation  exponentially. 

(3)  Battery  life  is  significantly  affected  by  the  ambient  tempera¬ 
ture,  the  too  low  EODV  or  too  high  EOCV,  and  the  excessive 
charge  and  discharge  rates.  Battery  life  shows  stronger  sensi¬ 
tivity  to  discharge  rate  and  EOCV  than  to  charge  rate  and  EODV. 
In  addition,  the  standard  charge  is  more  favorable  for  battery 
durability  than  the  fast  charge. 

This  degradation  model  alleviates  the  problem  of  poor  accuracy 
in  models  that  do  not  consider  coupling  and  demonstrates  a  prefer¬ 
able  adaptability  to  the  complexity  of  the  multi-factor  input.  The 
application  of  the  model  requires  only  the  calibration  of  coupling 
characteristics  and  the  input  of  real-time  operating  conditions, 
which  ensures  easy  applicability  to  batteries  on  electric  appliances 
and  vehicles. 
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